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ABSTRACT 


This thesis develops a multiple regression model using 
regional and unit characteristics to estimate commuting 
distances of U.S. Army Reservists. The data were obtained 
from a 1988 file established by the Defense Manpower Data 
Center containing locational and biodemographic information on 
238,174 enlisted reservists. A random sample of 91 reserve 
centers was selected for the analysis. 

The logistic and normal distributions were evaluated as 
possible candidates for fitting the commuting distance 
distribution. It was found that a power transformation of the 
fractional distance traveled fit both distributions quite 
well. Parameters for the two distributions are obtained 
through a method of maximum likelihood estimation. Finally, 
a multiple regression equation is used to estimate the 
parameters of the commute distance distribution as a function 
of reserve center and market characteristics. 

The results of the multiple regression equation provide 
the U.S. Army Recruiting Command with some important variables 


necessary to predict commuting distances. 
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I. INTRODUCTION 


A. GENERAL 

The Department of Defense's share of the national budget 
may soon be decreasing. As communism in Europe and Asia 
continues to crumble, funds may be reallocated from the Armed 
Forces to offset the national debt and bolster domestic and 
foreign expenditures. The U.S. Army is preparing for possible 
funding and personnel reductions. 

The Secretary of Defense has proposed four major reduc- 
tions in the 1991 Army budget that will save an estimated $3.3 


billion. These proposals are shown in Table 1 [Ref. l:p. 19}. 


TABLE 1 


PROPOSED REDUCTIONS IN THE ARMY'S 1991 BUDGET 


1. Deactivate 2 Army Divisions $1.2 


2. Cut Helicopter Improvement 
Program 







3. End Production of the M-1 
Tank 


4. Cancel Procurement of the 
Apache Helicopter 






Total: $3.3 














These sweeping proposals will shift a higher percentage of the 
Army's total defensive capability to the reserves. Currently, 
there are 18 active and 12 reserve divisions; 40 percent of 
the Army's fighting force is in the reserves. If the 
Secretary of Defense's budget proposal is approved, the force 
mix will be 16 active and 12 reserve by the end of 1991. The 
reserves will then comprise over 42 percent of the Army's 
defensive capability. Future possible reductions in 1992 and 
beyond foretell of increased reliance on the reserves. In 
short, the Army Reserve will continue to be an integral part 
of this country's defense forces. 

In light of the impending reductions, the reserves may 
expect to be affected in one of four ways. First, only 
active forces might be reduced, greatly increasing the current 
40 percent ratio of reserve to total Army units. Secondly, 
deactivated units might be transferred to the reserves, 
swelling the reserve strength and necessitating the opening of 
new reserve centers. Thirdly, the reserves will undergo troop 
reductions independent of active cutbacks. Lastly, the ratio 
of reserve units in the total Army force could remain at 
approximately 40 percent, implying that the reserves will 
suffer proportionally the same reductions as the active Army. 
While it is difficult to predict which optic, if any, will 
become reality, one can argue that the Army must be prepared 
either to establish new reserve centers or close down existing 


ones. 


Manning the reserve units is extremely important. A unit 


.cannot be expected to perform its mission unless it has 


trained, qualified personnel staffing critical positions. 
Units are considered undeployable if they cannot attain 
certain fill rates. 

The U.S. Army Recruiting Command (USAREC) is the proponent 
agency responsible for recruiting personnel to meet unit fill 
rates. An important objective of the USAREC is to identify 
market areas that will support new reserve units [Ref. 2:p. 
sae Most unit location studies have centered on supply 
variables such as civilian wages and unemployment (Ref. 2:p. 
11}. Little or no research has been done to establish a link 


between various supply variables and commuting distances. 


B. PURPOSE 

The purpose of this thesis is to analyze individual 
commuting distances of enlisted reservists and estimate a 
statistical relationship between commute behavior. and 
different supply variables. Such an estimated relationship 
can be used by the USAREC to analyze a geographical area and 
assess how it might support a new reserve unit. 

This thesis utilizes the logistic and normal distributions 
with suitable power transformations to model commuting 
patterns. The fitted distributions can then be used by the 
USAREC to determine the geographic extent of a recruiting 


market around a particular reserve center. 


C. BACKGROUND 

Commuting patterns have changed significantly over the 
past 100 years. In the period 1890 - 1920, people primarily 
commuted to the Central Business District (CBD). The main 
source of transport was the trolley [Ref. 3:p. 512}. The 
average commute distance using the trolley was 1.8 miles [Ref. 
4:p. 361). For convenience, houses were located within a 
short walking distance of the CBD or the trolley stop [Ref. 
3:p. 512]. 

After 1920, the introduction of the bus and private 


automobile increased commuting distances significantly (See 


Figure 1). 
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Figure 1 


Commute Distances: 1920-1985 











People had the freedom to disperse from transit lines, 
allowing cities to expand spatially. In the 1960's and 1970's 
an extensive highway and interstate system was developed to 
accommodate the dispersed population [Ref 4:p. 361]. 
Businesses began to move manufacturing and service sites out 
of the CBD to more suburban locations [Ref. 3:p. 515]. 
Because many businesses moved to the suburbs, commuting 
distances peaked at 11.1 miles in 1979 [Ref. 5:p. 3]. The 
trend for businesses to move from the CBD to the suburbs 
continued into the 1980's, as the average c.wuting distance 
decreased to 10.8 miles in 1985 [Ref. 6]. Therefore, one can 
infer that in the future commute distances of employees will 
stay the same or decrease slightly. 

Distances that personnel are willing to commute should be 
considered when assessing potential and existing Army Reserve 
Center locations. Research using commuting distances to 
predict manpower supply for the Naval Reserves was done by 
Beth Asch in 1985 [Ref. 7]. She contends that potential 
reservists will commute at most 100 miles. Therefore, she 
recommends that the Navy assess the geographical region within 
100 miles from a potential site to determine if the area can 
provide the necessary personnel to staff the ship(s). 
Likewise, the Army can use commuting distances to assess the 
feasibility of a geographical region to staff reserve units. 

A file containing location and biodemographic data on 


238,174 enlisted reservists es..cablished in 1988 by the Defense 





Manpower Data Center was used to develop this thesis. The 
individual commuting distances were estimated by using 
longitudinal and latitudinal data to obtain the Euclidean 
distance from the individual's home zip code to the unit's zip 
code. Therefore, the distances used in this thesis are 
estimates and not the actual commuting distances. Hence, one 
must assume that the relationship between the estimated | 
distance and the true distance remains the same within a 


reserve center all else being equal. 


D. ASSUMPTIONS 
Listed below are assumptions that were made prior to 
analyzing the data. 


1. Reserve centers are located at the centroid of the zip 
code region. 


2. Reservists reside uniformly throughout a zip code 
region. This is important because distances were 
computed by measuring the Euclidean distance from the 
centroid of the reserve certer zip code region to the 
centroid of the reservists' zip code. 


3. The average person will commute at most 100 miles. 
Therefore, commute distances in excess of 100 miles are 
not considered in this analysis. 


4. Commute distances of personnel residing in the same zip 
code region as the reserve center are greater than zero 
miles. This is important because the Euclidean 
distances in the database for these personnel equal zero 
miles. Therefore, through the use of a detailed road 
map, zip code boundz2zries were estimated and the distance 
across the zip code area was multiplied by 0.5 (0.5 was 
selected because this is the farthest possible commute 
distance to the reserve center). For example, if the 
distance across a zip code is 6 miles, then personnel 
residing in that zip code can travel at most 3.0 miles. 


5. Commute distances will vary from one region of the 
country to another. Therefore, 91 reserve centers were 


6 





selected using a random number generator from the 
Northeast, South, Midwest, and West. The number sampled 
from each region is proportional to the personnel in the 
reserves from that region. For example, since 
approximately 27 percent of reservists are located in 
the West, 25 of the 91 reserve centers (27 percent) are 
from the West. 
E. EXPLANATION OF DATA AND RESEARCH 

Appendix A contains Tables 2~5 listing the reserve centers 
selected in each geographical region. Some personnel did not 
have a home zip code, so they are not considered in the 
analysis. Also shown in each table are the number of 
personnel that commuted over 100 miles and the number living 
in the reserve center Zipcode. The necessary distance 
adjustments discussed in assumption number four are shown in 
Tables 6-9 in Appendix B. One should note that no more than 
10 percent of the reservists in a geographical region reside 
in the same zip code area as the reserve center. In addition, 
the analysis includes at least 91 percent of the distance data 
within 100 miles of each reserve renter. 

The next step is to evaluate the logistic and normal 
distributions and determine whether commuting distances can be 
accurately predicted. Because the data are not symmetric, a 
suitable power transformation is applied. Parameters for the 
logistic and normal distributions, in addition to the power 
transformation, are obtained through a method of maximum 
likelihood estimation using a numerical solver. The estimated 
parameters are substituted into the cumulative distribution 


functions to compute the appropriate percentile commute 
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distance. Finally, regression analysis techniques are 
investigated to determine if regional and unit characteristics 
can be used to predict the distribution parameters. These 
steps will be discussed in detail in the next three 


chapters. 


2) 


II. LOGISTIC DISTRIBUTION 


A. GENERAL 

Research characterizing part-time commute behavior for 
Army reservists was done by Laura D. Johnson and George W. 
Thomas using 30 reserve centers [Ref. 8]. Sampling from the 
1988 reserve data base, they conducted extensive exploratory 
data analysis to determine if there exists a probability 
distribution that adequately describes some power transfor- 
mation of the commuting distances. Relying primarily on the 
Kolmogorov-Smirnov statistic, they determined that the 
logistic distribution provided a good fit ([Ref. 8:p. 7]. 
Based on their results, the power logistic distribution was 
initially selected as a model of commuting distances of U.S. 
Army Reservists. 

The logistic distribution has previously been used as a 
growth curve and for demographic purposes [Ref. 9:p. 3}. It 
is symmetric and has a shape similar to the normal distribu- 
tion. The logistic has relatively "longer tails" meaning that 
it has greater variability about the mean [Ref. 9:p. 6]. For 
example, the inflection points on the standard normal curve 
are at +1 while the logistic distribution has its inflection 
points at +0.53 [Ref. 9:p. 6]. The relatively "longer tails" 
have their most significant impact on the fourth central 


moment. The value of the logistic is 4.2 vice 3.0 for the 








normal distribution [Ref. 9:p. 6]. Although the relatively 
"longer tails" have a considerable effect on the fourth 
central moment, there is a much smaller effect on the 
estimation of specific values of the cumulative distribution 
function (CDF). Thus, one could argue that it is appropriate 
to use either the logistic or normal distribution to predict 
the 75th and 90th percentile commute distances. Shown below 


is the logistic distribution [Ref. 9]. 


_ (yy ~@) 
_ e 
f(y) = wee) 2 
Bi1 +e B 
Ely] =a 
where ee, 
vary) = ce 


B. POWER TRANSFORMATION 

The commuting distances were examined to determine if a 
power transformation was necessary. One fast and efficient 
method is to plot the data on a quantile plot. As an example, 


the commuting distances for the reserve center at Long Beach, 


California are plotted in Figure 2. 
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Long Beach, CA N=151 


Logistic Quantile Plot Without a Transformation 


If the data do not fall close to the Y=X line, the data are 
not symmetric and will fail the Kolmogorov~-Smirnov statistic. 
(Ref. 10:p. 30] One way to get the data to lie close to the 
Y=X line is to apply a suitable power transformation to the 
commuting distances. 

The distance travelled by the ith member of a reserve unit 


is denoted by x, . The transformation used is [Ref. 8:p. 7]: 


il 





x,° if 6#0 
Ja } 
Inx, if 6=0 
If the random variable Y = u(X) defines a one-to-one corres- 
pondence between the values of X and Y, the equation y = u(x) 
will yield a unique value for x in terms of y. This relation- 
ship can be written as x = w(y). Therefore, the probability . 


distribution of Y is: 


g(y) = f£lw(y)] [J 


where |J| is the jacobian of the transformation. 
[Ref. ll:p. 143] 
To obtain the Jacobian, one takes the derivative of y with 


respect to x. 


Q@x.8-2 if 6#0 


ay . 
ax + iF 6=0 


Substituting into the equation for g(y), one obtains the 


probability density function (pdf) shown below. 


f£(x,;°) @x,5+ if 6#0 

Ge.) = 

: f(lnx,) + if 6=0 
a 


These equations will yield the commute distance for each 


individual after the power transformation. 


C. MAXIMUM LIKELIHOOD ESTIMATION 
A method of maximum likelihood was used to estimate the 
parameters of the logistic distribution. A numerical solver 


called the General Algebraic Modelling System (GAMS) was 


Py 


utilized to obtain values for a, §$, and 8 (call these values 
a, b, and t) [Ref. 12). Each reserve center has n 
| observations. The individual commute distances, x,, X,, ey 
x,, are used to obtain the following likelihood for each 
center (Ref. 13:p. 268]. 


Nn 


Ly(a,B) = [] 9(x,) 


2*1 


The purpose is to find estimates of a, b, and t that have the 


following desirable property [Ref. 13:p. 268]. 


: L. (a,b) 2 Ly(a,B) VW estimates of (a,f,8) 


These are called the maximum likelihood estimators (MLE) ana 
are obtained by taking the derivative with respect to a, b, 
and t and setting each equation equal to zero. Solutions to 
these equations were obtained numerically in GAMS. The 
derivations to obtain a, b, and t are shown in Appendix C. 
Obtaining optimal estimates can be time consuming because the 
bounds for a, b, and t must be established by trial and error. 
This process was done to each of the 91 reserve centers and is 


tabulated in Tables 10-13 in Appendix D. An example of one 


GAMS program and the output is in Appendix E. 











One way to easily verify the accuracy of the GAMS output 
is to plot the commute distances after applying the calculated 
power transformation (t). A quantile plot using the reserve 
center data from Long Beach, California is depicted in Figure 


3. 
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Figure 3 
Long Beach, CA N=151 


Logistic Quantile Plot With a Transformation 


The data lie very close to the Y=X line indicating that the 
fit is quite accurate. In addition, the transformed data 
easily pass Komogorov-Smirnov statistic, effectively 
validating the model. 

The a, b, and t values can be substituted into the 
logistic CDF to obtain estimated commute distances for the 


50th, 75th, and 90th percentile commute distances (See 
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Appendix F for calculations). The estimated commute distances 
from the logistic distribution versus the actual commute 


distances are in Tables 14-17 in Appendix G. 


D. RESULTS 

The power logistic distribution appears to be an accurate 
predictor of commuting distances. Estimated commute distances 
for the 50th percentile were very close to the actual 
distances. Let the prediction error be the actual distance 
minus the predicted distance. The absolute mean error is 1.8 
miles or less for each geographical region. Figure 4 is a box 
plot of the 50th percentile prediction errors. It is 
interesting to note that the greatest variability is in the 


West. Perhaps rural living accounts for this anomaly. 
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Figure 4 
50th Percentile Commute Distances 


Errors In The Logistic Predictions 
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The predicted commuting distances for the 75th percentile 
are shorter than the actual distances. This is not unexpected 
because of the "longer tails" discussed earlier. The 
prediction errors for the 75th percentile are depicted in 


Figure 5. 
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Figure 5 
75th Percentile Commute Distances 


Errors In The Logistic Predictions 


As expected, predictions in the Northeast are more accurate 
than the other regions. This is probability due to the highly 
developed, closely located cities in the Northeast. The 
‘spread of workers in the South and West and the resulting 
longer commute distances contribute to the greater variability 


in these regions. 
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Predicted commuting distances for the 90th percentile were 
very accurate in the Midwest and Northeast. Figure 6 
indicates that the West has much greater variability than the 


other three regions. 
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Figure 6 
90th Percentile Commute Distances 


Errors In The Logistic Predictions 


This is understandable because many communities in states such 
as Nevada and Colorado are long distances from the nearest 
reserve center. Several commuters traveling a long distance 
can skew predictions in the model. 

Overall, the logistic distribution is an accurate predic-~ 
tor of commuting distance percentiles. Rural living and the 
dispersion of communities might account for some of the 
greater variability in the South ane West. The more accurate 


Northeastern predictions may be a result of the highly 
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developed urban environment from Boston, Massachusetts to 
Washington, D.C. Another possible explanation for the 
smaller error in the Northeast is that most reserve centers 
have been purposely located to ensure that the 90th percentile 
commute distance is less than 50 miles. Predictions, 
particularly in the South and West, might have been effected 
by the elimination of distances in excess of 100 miles. — 
Lastly, it is conceivable that only 10% of reservists are 
willing to commute over 50 miles. One can surmise that if a 
center is not located within 50 miles of their residence, they 


seek alternate employment. 
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III. NORMAL DISTRIBUTION 


A. GENERAL 
As discussed in Chapter II, the logistic distribution is 
similar in shape to the normal distribution. Listed below are 
several characteristics common to both distributions (Ref. 
1l:p. 112]. 
1. It is continuous and symmetrical about the mean (p). 
2. The mode equals the mean. 


3. The curve has points of inflection at uw to (o is the 
standard deviation). 


4. The curve approaches the horizontal axis asymptotically 
to the left and right of un. 


5. The area under the curve is equal to l. 
Because of their similarities, one can speculate that commute 
distances under some suitable power transformation could also 
be modelled using the normal distribution. The hypothesis 
being tested is that the normal distribution will yield 
results that are as accurate as the logistic distribution. 


The normal distribution is depicted below [Ref. 13:pp. 90- 


92). 
-(vy-p)? 
f(v;)=—t-e 2 
V2no 
Elv;] =p 
Heese Var{v;] = 0° 
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B. POWER TRANSFORMATION 
The data were examined prior to applying a transformation. 
Again, a quantile plot was chosen because it adequately 


displays all of the data (See Figure 7). 
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Figure 7 
Long Beach, CA N=151 


Normal Quantile Plot Without a Transformation 


Because the data do not fit the Y=X line, an appropriate power 
transformation was applied. The transformation selected is 


Shown below [Ref. 14:p. 214]. 





OL 
if ne YT, 

V; = 0 
ln x if 6=0 


where x; is the distance travelled by the ith 


a 


member of the reserve unit. 
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This transformation is appropriate because it is continuous at 
O [Ref. 14:p. 215}. 

As discussed in Chapter II, if the random variable V = 
u(X) defines a one-to-one correspondence between the values of 
V and X, the equation v = u(x) will yield a unique value for 
xX in terms of v. This relationship can be written as x = 


w(v). Therefore, the probability distribution of V is: 
g(v) = flw(v)) |J| 


where |J| is the jacobian of the transformation. 
[Ref. ll:p. 143] 


To obtain the Jacobian, one takes the derivative of v with 


respect to x, . 


x, if #0 





3 - zt if 8=0 


Substituting into the equation for g(v), one obtains the 


following pdf: 





6 
f _ *) x,° if 8#0 


f (Inx,;) + if 60 
A: 


Following the same procedures as discussed in the logistic 


distribution, the likelihood function is obtained: 
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Iy(a,B) = [[ 9(x,) 


i=l 


a= p 


where 6 


iT 


6 = power transformation 


Substituting for g(x;) when 60 , one obtains: 





14(a,8) = TI te 20" (xi) 


isl onto 


After taking the natural logarithm of both sides one can 


simplify the equation to obtain: 





} >> (0-1) 1n (x;) 


i=l 


In Lg(a,B)=- Sinan - nino - 352 2. 5 (23 


Substituting for g(x,) when 6=0 , one obtains: 





‘i - (Qn xy-p)? =p)? 
Ly(a,B) = IT|  20F (=) 


One can take the natural logarithm of both sides and obtain: 





Nn n 
In Lela, B) =-5 In 2x ~nln - (In x;-p)? - ln x, 
i=l 


i=} 


After obtaining the estimated parameters from GAMS, these 
equations will yield the commute distance distribution for 


each reserve center. 
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GAMS was used to implement a method of maximum likelihood 
to estimate the parameters of the normal distribution. The 
estimators for a, £, and ® are labelled m, s, and t. An 
example of one GAMS program and the associated output is in 
Appendix H. Values for each of the 91 reserve centers were 
computed in a Similar manner and are tabulated in Tables 18-21 
in Appendix I. 

The GAMS model to obtain the estimates can be verified by 
graphing the data. After applying the power transformation t, 


the data were graphed on a quantile plot (Figure 8). 
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Figure 8 
Long Beach, CA N=151 


Normal Quantile Plot With a Transformation 
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Clearly, the transformed data lie very close to the Y=X line. 
In addition, the Kolmogorov-Smirnov statistic is satisfied 
providing further proof that the transformation is accurate. 

The values of m, s, and t are substituted into the normal 
CDF to obtain estimated commute distances for the 50th, 75th, 
and 90th percentiles. Calculations are in Appendix J. The 


estimated commute distances are in Tables 22-25 in Appendix K. 


C. RESULTS 

The normal distribution with the appropriate power trans- 
formation is an excellent predictor of commuting distances. 
Figure 9 is a box plot of the 50th percentile prediction 


errors. 
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Figure 9 
50th Percentile Commute Distances 


Errors In The Logistic (L) And Normal (N) 
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The largest absolute mean error is 2.1 miles in the South. 
Prediction variability is greatest in the South, a phenomenon 
that was also observed with the logistic distribution. 
Estimators are skewed right in the Northeast and South but are 
skewed left in the West. One possible explanation for 
predicted distances exceeding actual distances in the west is 
the geographical dispersion of communities. Many western 
communities are close together because of the terrain and 
water availability. Therefore, several miles may separate 
Clusters of communities. Consequently, the commuting flow 
consists of many different groups of personnel travelling 
almost the same distance. 

Predicted commuting distances for the 75th percentile are 


depicted in Figure 10. 
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Figure 10 
75th Percentile Commute Distances 


Errors In The Logistic (L) And Normal (N) 
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As expected, prediction errors in the Northeast are less than 
other regions. Although variability is greater in the West, 
the mean absolute error is only 4.7 miles. This relatively 
small error indicates that the normal distribution provides a 
good estimate. The mean absolute error is 5.2 miles in the 
Midwest. Even though this is higher than the error in the 
West, the variability is much less, implying that the normal 
Gistribution will provide relatively consistent predictions. 

Variability and skewness in the 90th percentile 
predictions for each region are similar to predictions for the 
75th percentile. Figure 11 shows that predictions are quite 


accurate in the Northeast and Midwest. 
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Figure 11 
90th Percentile Commute Distances 


Errors In The Logistic (L) And Normal (N) 
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Rural living and the associated spatial movement away from 
large metropolitan areas contributes to the greater 


variability in the South and West. 


D. COMPARISON OF THE NORMAL AND LOGISTIC DISTRIBUTIONS 

The normal and logistic distributions with suitable 
transformations have been used to predict commute distances. 
Figure 12 depicts prediction errors for the 50th, 75th, and 


90th percentile commute distances for 91 reserve centers. 
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Figure 12 
Percentile Commute Distances 


Errors In The Logistic (L) And Normal (N) 


The 50th percentile predictions are nearly identical for 
both distributions. A slight difference in variability 
between the distributions is evident in the 75th percentile. 
Both distributions are skewed right, with the normal yielding 
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more consistent predictions. The 90th percentile predictions 
for both distributions are nearly identical. It appears that 
the "longer tails" of the logistic distribution have little 
effect when compared to the normal distribution. Because the 
predictions from both distributions are nearly identical, one 
can select either distribution to estimate commuting 
distances. 

Although either distribution can adequately forecast 
commuting distances, the normal distribution is preferred for 
several reasons. The MLE's for the normal distribution are 
faster to compute in GAMS than the MLE's for the logistic 
distribution. Secondly, the MLE's for the normal distribution 
are asymptotically normally distributed. [Ref. 13:p. 272). 
In order to verify that the MLE's possess normal distribution 
properties, quantile plots with Kolmogorov-Smirnov bounds were 
evaluated for m, s, and t. Shown below in Figure 13 is the 


quantile plot of the t values. 
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Figure 13 
Quantile Plot of t Values 
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The values clearly lie along the Y = X line and are well 
inside the Kolmogorov-Smirnov bounds, indicating that the t 
values are normally distributed. Similar quantile plots were 
evaluated with identical results for the m and s values. 
Thus, the MLE's are consistent with the normal distribution. 
Conversely, the MLE's for the logistic distribution may not be 
asymptotically efficient because the transformation is 
discontinuous when @ = 0. In summary, because the computed 
power normal distribution is continuous and faster to compute 


it is preferred over the power logistic distribution. 
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IV. MULTIPLE REGRESSION 

A. GENERAL 

The objective of this chapter is to establish a multiple 
regression model using local geographical and unit character- 
istics that will best predict the values of u, o, and 6, based . 
on the estimators m, s, and t found in Chapter III. Listed 
below are geographical and unit characteristics that were 
tested as explanatory variables in the regression models. 

1. Geographical Characteristics 


a. The number of males in category III-A' and above from 
17-29 years old residing in the reserve center county. 


b. The median household income in the reserve center county 
Givided by the average annual military wage at the 
reserve center. 


c. Unemployment rate in the reserve center county. 


d. The percentage of full time civilian workers that 
commute to work outside of their residential county. 
This variable is a surrogate measurement for the 
propensity to commute. 


2. Unit Level Characteristics 


a. The percentage of male personnel stationed at the 
reserve center. 


b. The number of prior service personnel (both men and 
women) serving at the reserve center. 


‘potential recruits that are classified in category III-A or 
higher must meet the following educational and _ testing 
requirements: High School Degree (HSDG) or Currently in High 
School past the 10th grade (CIHS) or Alternative Credential Holder 
(ACH) and scores of at least 50 on qualifying Armed Services 
Vocational Aptitude Battery Tests (ASVAB) [Ref. 15 :p.14]. 
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In addition to the geographical characteristics, the model was 
tested by partitioning the United States into nine census 
regions. The census regions and the states assigned to each 
region are listed in Appendix L. The model was also tested by 
unit type. Each reserve center was categorized as combat, 
combat support, combat service support, or medical based on 
the largest type of unit assigned to the reserve center. A 
listing of the different types of units and which of the four 


categories they were assigned to is in Appendix M. 


B. ASSUMPTIONS 

Prior to developing and testing the multiple regression 
model, several assumptions were made concerning the geograph- 
ical and unit characteristics. These assumptions are listed 


below. 


1. Salaries increased four percent annually during the 
period 1979 to 1988. This is important because the most 
current tabulated figures for median household income by 
county is 1979 [Ref. 16]. In order to arrive at an 
income for 1988, the 1979 income figures were adjusted 
four percent annually. 


2. Propensity to commute can be captured by using the 
percent of workers that commute to work outside of their 
residential county. The data for each county were 
obtained from the Bureau of the Census [Ref. 17]. 


3. The potential recruiting population for the reserves is 
adequately represented by the number of male category 
III-A's residing in the reserve center county. 


4. Recruiting priorities for a reserve center are oriented 
toward filling the largest unit. Therefore, a reserve 
center can be categorized into one of four unit types 
(1.e., combat, combat support, combat service support, 
and medical) if fifty-one percent or more of the 
assigned personnel belong to that type of unit. 
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5. Regional differences in lifestyles are captured by 
dividing the county into nine census regions. 


6. The percentage of civilian workers commuting outside of 
their residential county remains unchanged during the 
period 1986 to 1988. This is important because the most 
current data on the percentage of workers commuting 
outside of their residential county is 1986 [Ref. 17]. 

C. THE GENERAL MODEL 
Multiple linear regression is used to describe a dependent 


variable by several independent explanatory variables. The 


general model is shown below [Ref. 10:p. 245]: 


Ye byt Udy xy + e, for iz=iton 
#1 


where y, are the n observed values of the response variable, 
X,, are the n values of the kth explanatory variable, 
e, are the normally distributed error terms, 
and b,, b,, ..., b, are the unknown regression 


coefficients. 
In order to use multiple regression the following assumptions 
must be satisfied [Ref. 18:p. 412}. 


1. The relationship between the dependent and independent 
variables is linear. 


2. The x values are observed without error. In addition, 
X values are known constants associated with the random 
Variable y. 

3. The y values are mutually independent random variables. 


4. The variance of y is constant for all x values. 








5. 


D. 


e. is normally distributed with a zero mean and 
constant variance. 


THESIS MODEL 


The model for this thesis is shown below. 


yp = b+ Edy Ky e; for 1 = 1 to 91 


where x, = male category III-A population in the county 

X, = percentage of males assigned to the reserve center 

X, = annual median household income divided by the 
average annual military income at the reserve 
center | 

x, = unemployment rate in the reserve center county 

x, = total number of male and female prior service 
personnel divided by the total numher of personnel 
assigned to the reserve center 

XxX, = percentage of workers that commute to work outside 


of their residential county 


XxX, = the census region or the unit type. 


The model assumptions are listed below. 


1. 


The relationship between pp, o, and @ and the selected 
independent variables is linear. 


Values for each of the independent variables are known. 


The values of m, s, and t are mutually independent 
random variables. 


The variances of m, s, and t are constant for the 
seven explanatory variables. 
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5. The distributions of m, s, and t are normal. 
The implicit null hy*othesis is that the variable coefficients 
equal OQ. 

The model was run using the General Linear Model Procedure 
(Proc Glm) in SAS [Ref. 19]. Initially, it was tested with 
the first six explanatory variables. The seventh explanatory 
variable, census region or unit type, was not considered in 
order to determine if the model can be applied nationwide 
without regard to census region or unit type. Statistics from 
the analysis of variance tables for m, s, and t are summarized 
in Table 26 in Appendix N. To insure that the assumption of 
normality was valid, the residuals for m, s, and t were 
evaluated using quantile plots. Shown below in Figure 14 is 
the quantile plot with Kolmogorov-Smirnov bounds of the 


residuals for the t values. 
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Figure 14 


Plot of Theta Residuals 


The residuals lie along the Y=X line and clearly fall inside 
the Kolmogorov-Smirnov bounds. Therefore, the assumption of 
normality is valid when the dependent variable is t. Similar 
quantile plots were evaluated for the m and s residuals but 
these failed the Kolmogorov-Smirnov test. ‘he reasons for 
failure are the m and s estimates from the reserve center at 
Fond du Lac, Wisconsin. Referring to Table 24 in Appendix K, 
one observes a great disparity between the actual and 
estimated commuting distances for the 75th percentile. 
According to Table 4 in Appendix A, 88 out of 160 distance 
data points were adjusted because the reservists reside in the 


Same zip code region as the reserve center. This inordinately 
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high number of adjustments is the most likely cause for the 
Tne Mend Ss residuals, 
Because over 55 percent of the data points were modified, it 
is not unreasonable to drop the m and s values. After 
elimination, the residual plots pass the Kolmogorov-Smirnov 
test. Therefore, one can conclude that the assumption of 
normality is also valid when the dependent variables are m and 
s. Model results for m, s, and t are shown in Tables 27-29 in 
Appendix O. 

The model was then tested including the nine census 
regions as the seventh explanatory variable. Census regions 
were selected as a test variable rather than the four 
geographical regions (i.e., northeast, midwest, south, and 
west) because census regions are commonly used to describe the 
United States. Model results for m, s, and t are shown in 
Tables 30-32 in Appendix O. 

Final model testing included unit type as the seventh 
explanatory variable. Model results for m, s, and t are shown 


in Tables 33-35 in Appendix P. 


E. MODEL RESULTS 

1. The initial model pinpoints the variables that USAREC 
should carefully evaluate when assessing a geographical area. 
When applied nationwide, the full model R-Squared values for 
m, s, and t are between 0.1104 and 0.1619 (Table 26). These 
are not particularly high values, indicating that the model 
does not predict a significant portion of the dependent 
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variables' variance. For example, the model predicts at most 
16.2% of the variance in the t values. One possible reason 
for the low R-Squared values is the fact that much of the data 
is estimated. For example, the commute data consists of 
Euclidean distances and not the actual distances reservists 
travel to their reserve centers. Another possible reason is 
that many important factors that impact on an individual's 
decision to commute may not have been captured in the data. 
For example, the road network in an area or the cost of fuel 
has not been captured in this model. However, upon examining 
Tables 27~29 one can identify those independent variables that 
have the greatest influence when predicting values of m, s, 
and t. One important statistic when evaluating the model is 
the last column (the probability that the variable is greater 
than the absolute value of the t statistic). Low probabil- 
ities indicate that the independent variable has a significant 
impact on the dependent variable. For example, in Table 27 
in Appendix N the variable, propensity to commute, equals 
0.0097. If our significance level is 0.05 (the value of 
a=0.05) then propensity to commute is a significant variable. 
Conversely, the p-value for unemployment is greater than 
0.05. Hence, unemployment is not a significant variable. 
Therefore, the significant variables are those that have a 
probability less than 0.05. 

The parameter estimates in Tables 27-29 in Appendix O can 


be used to construct the multiple regression equations. For 
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example, the regression equation for the dependent variable t 


trom Table 27 is snown beiow. 
t = -0.7518174 + 0.6258805 x, + 0.0063553 X, 


Equations for m and s can be constructed in a similar manner 
from Tables 28 and 29. 

2. The second series of multiple regression equations . 
were calculated using the nine census: regions as a dummy 
variable. The analysis was conducted in SAS using Proc GLM 
commands with census regions as a class variable [Ref. 19:p. 
434-506]. Using the value in the last column (the probability 
that the variable is greater than the absolute value of the t 
statistic), the mountain census region is a significant 
predicator for m and t when a equals 0.05. The variable s is 
best predicted by the mountain census region but it is not 
Significant unless a is greater than 0.0532. Except for the 
mountain census region, it appears that there is no difference 
in s for nine census regions. 

3. Unit type was used as a dummy variable in establishing 
the last three multiple regression equations. Again, the t 
statistic was used to determine if unit type was a significant 
variable. The lowest value was 0.4191 for variable "combat 
support" located in Table 33 in Appendix Q. Because all the 
t statistic values exceed 0.05, unit type, as measured in this 
model, does not appear to be not a significant variable for 


predicting values for m, s, and t. 
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F. CONCLUSIONS 

The purpose of this thesis was to analyze individual 
commuting distances of enlisted U.S. Army Reservists and 
estimate a statistical relationship between commute behavior 
and different market variables. After analyzing commute 
distances from a sample of 91 reserve centers, it was 
determined that either the logistic or normal distribution 
with a suitable power transformation could accurately predict 
commuting distances. Various combinations of regional and 
unit characteristics were tested with the goal of developing 
a multiple regression equation to predict the parameters of 
the normal distribution. Although the R-Squared values were 
low, the model does identify civilian propensity to commute 
and percentage of male personnel in the unit as the two most 
influential independent variables. 

Further research should focus on several independent 
variables not available in the 1988 data base established by 
the Defense Manpower Data Center. Four variables that would 
be invaluable are listed below. 

1. Cost of fuel in the various census regions. 

2. Individual commuting times. 

3. Local road conditions. 

4. Road access to each reserve center. For example, a 
reserve center located adjacent to an interstate might 


induce individuals to commute farther than one located 
off a secondary road. 
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Once these and other pertinent variables are included in the 
regression model the ability to assess a geographical region 


may be greatly enhanced. 
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APPENDIX A 
TABLE 2 


NORTHEASTERN UNITED STATES 


[aon momen [orm] [ee [a [ae 
Mi. Zip 
East Windsor, CT 06088 | s30 | 23 | 19 | 6 | 
West Hartfora, cr_| os110 | 575 | 28 | 25 | 5 _ 
[Boston mA | ozz10_| sea | 38 | 52_ | a2 
[Lawrence MA__—|_o1ea3,_ | asa | a4 | 7 | 30 
02131 1 
04402 
04930 
08067 
14424 
13617 

104 

04 a 

55 

119 

19 

50 


ho 
~~ 
— 
~ Oo Ww jx] 
_ 
Ce) 


i) 
- 


2 


N) 
OF 


21 


N 


Elizabethtown, NY 12932 104 
Glens Falls, NY 12801 
Allison, PA 15601 | 304 


3 
3 






4 
2 
3 
3 
9 


ra) 
~J 


17 


16504 | 355 | a 
Huntingdon, PA 16652 119 _ 
Philadelphia, PA 19154 1 


15401 3 
5368 267 


_ 
NO 


7 
3 
4 


~J] 


Uniontown, PA 


Le] 
oO 
or 
~~] 


N 
NM 
~ 
WwW 
N 
> 





Percent of data used: 91 
Percent unknown: 05 
Percent greater than 100 miles: 04 


Percent in the same zipcode: 06 
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TABLE 3 


SOUTHERN UNITED STATES 


Center Location | Zipcode Missing | >100 Mi. | Same Zip 
35601 2 18 
3 
| 3s | to 
10 
12 


AL 
Dothan, AL 36302 a 
AR 


72702 
72730 
72204 
72301 
19973 


Tallahassee, FL 32304 365 11 
Bardstown, KY 40004 121 
Slidell, LA 70459 
Hattiesburg, MS 39401 


paz. 
39209 | 509_ 
27409 | 290_ 
28401 


1 
1 
1 


O* 
_ 


te 
~] 


56 
17 
87 
120 
183 


_ 
OV 
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© 
oO 





~~ 
On 


te 
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We) 
NO 
oO | & 
Ht 10 


~ 
8 


74003 
Greenwood, SC | 29646 


Carrollton, TX 75006 
Ft. Bliss, TX 79906 
Houston, Tx | _77054_| 1080 





> 
LO! 


- Oo 
WwW | © a 


~ 
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TABLE 3 (CONTINUED) 





center Location 
a7 | oe | 
115 
Ft. Story, VA 681 
U 
VA 







Richnond, VA | 23220 
24153 
charleston, WW | 25313 







Percent of data used: 91 


Percent unknown: 03 
Percent greater than 100 miles: 06 


Percent in the same zipcode: 10 
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TABLE 4 


MIDWESTERN UNITED STATES 


Center Location | zipcode | n_| Missing | >100 Mi. 
s0010_| 158 23 


po 
7 
ss 
144 Sosa 





















196 

Chicago, IL | so629 | sas] a4 | a7 

[Granite city, 11 | 6zos0 |sea| 6 | a7 

quincy, m | 6230 far] an | |e 
‘springfield, 11 | 62703 | 102] a 

Fort wayne, Iv | 46809 | 522] 8 | a9 | a3 
cary, IN| 46404 | 208 

307 

Scottsburg, IN 136 ae 

306 


Grand Rapids, MI 






Washington, MO 63640 1 
columbus, OH agaa5_ [7a] a2 | 48 | 0 
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TABLE 4 (CONTINUED) 


>100 Mi 
are 
ror | oo 


Wausau, wr | _54402_| 225 | 3 
roman | oz |e] 


Percent of data used: 91 





Percent unknown: 03 
Percent greater than 100 miles: 06 


Percent in the same zipcode: 10 


45 





TABLE 5 


WESTERN UNITED STATES 


__center Location | zipcode | _n_| missing | >100 mi. | sane Zip | 
Tucson, AZ [esis | sez | 49 | 93 | 29 
Bakersfield, ca__| 93301 | 191 | 22 | 20 | 5s 
[chico, ca 95926 | 224 
concord, CA _sasi9 | ago | 2 | 2 | 4 
Dublin, CA 
El Monte, CA | 91733 | 639 | 4 14 


90822 





Long Beach, CA 


Redding, CA 96003 a2 | 
San Pablo, CA 94806 431 


Wd 
© 


te 


NO a 
Mire IN fr IN Jw 


1 
1 
Upland, CA 91786 1 


39 
184 
112 
Vallejo, CA 94589 
23 
89 
137 
124 
188 


6 

Van Nuys, CA 91403 

Aurora, CO 80011 
3 


3 
7 


> 
td 


of 


15 
2 


Aurora, CO 80045 
Denver, CO 80225 | 389 | 


1. 
6 
0 
0 
2 
7 
0 


aad 


ois > 
io | > he ifs] 7h 


Hayden Lake, I 83835 


D | 10 
Great Falls, MT 59403 Fg 
Helena, MT 59601 ra 


Las Cruces, NM 88001 piss} 3 17 


14 


56 
oboe 
nem 
| 20 
43 

1 

1 

1 


3 
0 
3 
3 
5 
0 
4 


TABLE 5 (CONTINUED) 





Percent of data used: 91 


Percent unknown: 03 
Percent greater than 100 miles: 06 


Percent in the same zipcode: 10 
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APPENDIX B 
TABLE 6 


NORTHEASTERN UNITED STATES 


jeest Windsor, of 2:80 
est mrton,ct [nas 
[Rosiingale, wa | 2.00 



















4.50 
3.00 
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TABLE 7 


SOUTHERN UNITED STATES 





Center Location 
East Camden, 
El Dorado, 
Little Rock, AR 
West Memphis, AR 
I 
S 


L 
AL 
AR 
AR 
| 
NC 
OK 


Wilmington, 





Greenwood, SC 





Carrollton, TX 


Ft. Bliss, TX 
Houston, TX 
















Mileage Adjustment 


6 


oi 
Lo) 


6.50 
6.00 
5.00 


6.50 

3.00 

3.00 
12.50 





49 


TABLE 7 (CONTINUED) 


Center Location Mileage Adjustment 





ae... ae 
[eictoria, wx | 5.00 
Tre, story, a | s.00 
[eictona, va | _#.00 
|. 

= 


T 
Charleston, 
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TABLE 8 


MIDWESTERN UNITED STATES 


Center Location 
Council Bluffs, 6. 
I 





Mileage Adjustment 





Oo 
© 


IA 
Decorah, IA 
5.00 
6.00 
Granite City, IL 4.50 

IN 

I 

N 


7.50 
5.00 
2.50 
2.00 
5.00 
4.50 
3.00 
8.50 


5.00 
6.50 
5.00 


A 
f 
IA 
IL 
Quincy, IL 
i 
IN 
N 
IN 
MI 


Columbia, MO 





Washington, MO 
Columbus, OH 





51. 


TABLE 8 (CONTINUED) 


Center Location Mileage Adjustment 
Toledo, OH 


‘wmusen, wr | 6.00 





52. 





TABLE 9 


WESTERN UNITED STATES 


2.00 
3.00 
Aurora, CO 2:00 


Aurora, 


6.00 
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TABLE 9 (CONTINUED) 


Mileage Adjustment 





Center Location 











5.00 
fogden, ur] 
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APPENDIX C 


In the GAMS program for the logistic function, the equation 


, LIKE was used if 8#0 and LIKE2 was used if @=0 . Listed below 


is the derivation of each equation. 


Equation Like 


Ly(a,B) = [] 9(x;) 


isi 


Substituting in the equation for g(x,;) , one obtains: 


. (x ,°~«) 
Ig(a,B) = a 6 x? 


a4 _ (xsP-a) \* 
B\l +e 


The equation can be manipulated into the following forn: 
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By taking the natural logarithm, one obtains: 


(S> xi- na) n 
Ly(a,B) = enaaee ‘amir + (nin 6) + (8-1) }> x, 
im1 
n __(xg-a) 
-ninB-2y)lIn(i+e 6 ) 
i=1 


Equation Like2 


_ (in x,P-a) 
e 


Ly(a,B) = [J 


ae1 


o 
5 4 << (In x, =| xX; 


n 


(¥7 in xj’ na) 
e B 1 
n _ (ln x/%-a) \2 
6%) [use II *: 


1*1 





Ly(a,B) = 
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Taking the natural logarithn: 


y $i 
(S" in x,°- na) _ (in x,f-a) 


Iy(a,p) = -—£_——- n In B- 2¥ ink +e B } - Sing, 


p i=i ixl 


GAMS numerically solves equations Like and Like2 for a, £6, and 8 


(these estimates are referred to as a, b, and t in Chapter II). 
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APPENDIX D 
TABLE 10 


WESTERN UNIT_D STATES: LOGISTIC DISTRIBUTION 


-~0.3203 


Center Location Zipcode 
85713 
93301 


0.335 
0.162 


0.082 
0.023 
0.832 0.036 
1.285 0.083 
0.610 | 0.080 


Tucson, AZ 





1.422 
1.240 
0.540 


Bakersfield, CA 





Chico, CA 


Concord, CA 94519 
Dublin, CA 94566 


El Monte, CA 91733 
Long Beach, CA 90822 
Redding, CA 


91786 
94589 
91403 
80021 
80045 


Hayden Lake, ID 








0.631 





0.748 0.053 
1.020 
0.837 0.033 
1.606 0.184 
1.148 
1.444 0.125 


1.824 ~0.3918 
3.241 0.621 -0.4692 
0.507 ~0.4229 


0.496 | -0.4481 
0.998 -0.5384 


0.430 ~0.4140 
0.967 ~0.5649 
1.018 0.005 ~0.0088 
1.502 0.124 -0.1740 
2.333 0.452 ~0.3858 


83835 
59403 
59601 





Great Falls, MT 





Helena, MT 
Las Cruces, NM 88001 
Reno, NV 89502 
Eugene, OR 97402 


84321 
84407 
98404 
99216 





2.479 
4.329 
2.400 








O1 
© 


TABLE 11 


MIDWESTERN UNITED STATES: LOGISTIC DISTRIBUTION 


Center Location 


Council Bluffs, IA 





Decorah, IA 


Dubuque, IA 
Chicago, IL | 


Granite City, IL 


Springfield, IL 
Fort Wayne, IN 


IN 
Lake Station, IN 
Bay City, MI 
Grand Rapids, MI 


Duluth, MN 
Fergus Falls, MN 


Columbia, MO 
Columbus, OH 
















Washington, MO 
Toledo, OH 
Beloit, WI 


Fond du Lac, WI 


Wausau, WI 





Zipcode 
50010 
51501 
52101 
52001 
60629 
62040 
62301 
62703 
46809 
46404 
46405 
47170 
48706 
49503 
55802 
56537 
65201 
43215 
63640 
43606 
53511 
54935 
54401 


59 


a ae 
1.234 
1.587 | 9.134 
0.800 | 0.054 


0.954 0.014 
1.652 0.214 


3.995 0.929 
~2.670 0.500 


1.034 0.011 


t 
-0.0452 
-0.4813 
0.2317 
-0.4219 
-0.2030 
-0.5208 
0.2152 
-0.0989 
-0.1735 
-0.1713 
-0.2453 
0.4770 
0.1328 
0.1122 
-0.1019 
0.5511 
-0.3507 
-0.1564 
0.0302 
~0.1933 
-0.6248 
0.0000 
-0.0199 


TABLE 12 


SOUTHERN UNITED STATES: LOGISTIC DISTRIBUTION 


El Dorado, AR 





Ft. Story, VA 


Richmond, VA 


Salem, VA 





Charleston, WV 


| zipcode | a | ob | 

0. 6003 
0.1465 
0.1000 
0.1075 
-0.2511 
0.496 
| 19973 | 0.677 | 0.062 


-0.5349 
0.2633 
-0.0015 


| 0.496 | 
0.3335 
rosso | 0.406 | 0.092 | 0.7523 
0.0152 
0.5219 
27409 0.970 0.0192 
| 9.480 _| 


28401 2.675 0.480 ~0.4495 


-0.1179 
0.2228 
-0.4593 
-0.0729 
0.1928 
-0.1250 
-0.3200 
~0.1647 


SS ET Sa a TT ET ERI EF 


74003 0.120 | 0.6938 


OV 
om) 





TABLE 13 


NORTHEASTERN UNITED STATES: LOGISTIC DISTRIBUTION 


center Location | zipcoee | a | ob | ot 
East Windsor, CT 06088 0.097 | -0.1548 
Middletown, CT 06457 | 0.962 _ 








_ 0.014 _| 0.0225 

West Hartford, cr__| os110__| 0.927 _| 0.022 | _0.0417_ 

ftawrence, MA |__ 02843 | 0.602 | 0.069 | 0.2394 

Roslindale, mA | 02232 |_4.912_| 1.005 _| -0.4952_ 

1.626 

[Pedrickton, nJ__ | __ogos7_|_0.681__| 0.059 | 0.2229 
15601 





15601 | 0.602_| 0.072 | 0.2267_ 
|0.020_| -0.0542_ 
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APPENDIX E 

A GAMS program was used to numerically solve the MLE's (a, b, 
and t) using the logistic distribution for every reserve center. 
The input was formatted by using a fortran program. An example of 
the GAMS program is shown below. The reserve center depicted is 
located at Bakersfield, California. 

PROGRAM 

SOFFSYMLIST OFFSYMXREF 
OPTIONS SOLPRINT = On, LIMCOL = 0, LIMROW = 0; 
sets 


Sinclude file set; (this file assigns a number to each commuter; 
in this example, the numbers are 1 to 149) 


scalar count 


Sinclude file count; (this file gives a scalar for the number of 
commuters to the reserve center; in this 
example, count equals 149) 


parameter X(I) 


Sinclude file out; (this file is a listing of the commute 
distances) 


variable 
value (this is the value of the maximum likelihood function) 
al (this is alpha) 
be (this is beta) 
th; (this is theta, the value of the power transformation) 


positive variable be; 


be.lo = 0.010; (the range for beta is: 0.00 < beta < 1.0) 

be.up = 1.000; 

al.lo = 0.010; (the range for alpha is: 0.00 < alpha < 5.0) 
al.up = 2.000; 

th.lo = 0.010; (the range for theta is: 0.00 < theta < 1.0) 
th.up = 1.000; 

al.l = 0.180; (these four lines establish a starting point for 
be.l = 0.280; the algorithm; the "1" means level) 

th.l = 0.250; 

value = 0.000; 
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equations 
LIKE2 
LIKE1; 


LIKE1.. value =e= 
count*(al/be) - SUM(I, (X(I) **th) )/be 
- count*LOG(be) - 2*SUM(I,LOG(1 + EXP(-((X(I) **th) ~- al)/be))) 
+ count*(LOG(ABS(th)) + (th-1) *SUM(I,LOG(X(I))); 


Model Stat /LIKE/; 

Solve Stat using DNLP maximizing value; (this directs the solver 
to find the maximum value for the function 
LIKE using the Discontinuous Non-Linear 
Program algorithm) 


Display be.l, al.l, th.1, value.1; (this displays the values for 
alpha, beta, theta, and the function LIKE) 


be.lo = 0.006; (the next series of lines sets new bounds and 
be.up = 1.000; starting values for alpha, beta, theta and value) 
al.lo = 0.010; 

al.up = 3.600; 

th.lo = -1.000; (the range for theta is: -1.00 < theta < 0.00) 
th.eup = -0.001; 

al.l = 0.250; 

be.l = 0.180; 

th.l = -0.200; 


Solve Stat using DNLP maximizing value; 

Display be.l, al.1l, tn.1, value.]; 

LIKE2.. value =e= 
- count*~L0OG(be) - SUM(I,LOG(X(I)))/be + count*(al/be) 
- SUM(I,LOG(X(I))) 
~ 2*SUM(I,LOG(1 + EXP(-(LOG(X(I))-al)/be))); 


Model Stat2 /LIKE2/; 


be.lo = 0.250; (the next series of lines sets the bounds and 
be.up = 1.000; starting values for alpha, beta, and theta for 
al.lo = -4.000; the equation LIKE2) 

al.up = -0.010; 

be.l = 0.280; 

al.l = -0.280; 

value = 0.000; 
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Solve Stat2 using NLP maximizing value; (this directs the solver to 
use the Non-linear Programming algorithm) 


Display be.1, al.l, value.1l; (the value of theta is zero so it 
is not being displayed) 


be.lo = 0.2500; (the next series of lines sets new bounds and 
be.up = 2.0000; starting values for alpha and beta) 

al.lo = 0.0005; 

al.up = 2.5000; 

be.l = 0.2800; 

al.l = 2.0000; 

value.l = 0.0000; 


Solve Stat2 using NLP maximizing value; 


Display be.1, al.1l, value.1l; 


Output 

LOWER LEVEL UPPER MARGINAL 
EQU LIKE1 320.2040 320.2040 320.2040 1.0000 

LOWER LEVEL UPPER MARGINAL 
VAR Value -~INF 66.8657 +INF ‘ 
VAR al 0.0100 0.9710 2.0000 EPS 
VAR be 0.0100 0.0100 1.0000 -34.2779 
VAR th 0.0100 0.0133 1.0000 EPS 
REPORT SUMMARY: 0 NONOPT 


0 INFEASIBLE 
O UNBOUNDED 


0 ERRORS 
VARIABLE be.1i = 0.010 
VARIABLE al.l = 0.971 
VARIABLE th.l = 0.013 
VARIABLE value.l = 66.866 
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LOWER 


EQU LIKE1 320.2040 


LOWER 
VAR Value -INF 
VAR al 0.0100 
VAR be 0.0060 
VAR th -1.0000 


REPORT SUMMARY: 


VARIABLE be.1] 
VARIABLE al.1l 
VARIABLE th.1l 


VARIABLE value.1l 


LOWER 


EQU LIKE2 320.2040 


LOWER 
VAR value ~INF 
VAR al -4.0000 
VAR be 0.2500 


REPORT SUMMARY: 


VARIABLE be.] 
VARIABLE al.l 


VARIABLE value.1l 


LEVEL 
320.2040 

LEVEL 
68.9552 

1.4218 


0.1623 
-0.1535 


0 NONOPT 
O INFEASIBLE 
0 UNBOUNDED 
0 ERRORS 
0.162 
1.422 
“0.153 


68.955 


LEVEL 


LEVEL 


67.1863 
“2.2136 
0.7689 


0 NONOPT 
0 INFEASIBLE 
O UNBOUNDED 
0 ERRORS 
0.769 
“2.214 


67.186 
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320.2040 


UPPER 
320.2040 
UPPER 
+INF 
3.6000 


1.0000 
-~0.0100 


UPPER 


320.2040 


UPPER 


+INF 
-~0.0100 
1.0000 


MARGINAL 
1.0000 
MARGINAL 
EPS 
EPS 
EPS 


MARGINAL 
1.0000 
MARGINAL 


EPS 
EPS 


LOWER 


EQU LIKE2 320.2040 
LOWER 
VAR Value ~INF 
VAR al 0.0005 
VAR be 0.2500 


REPORT SUMMARY: 


VARIABLE be.1l 


VARIABLE al.l 


VARIABLE value.1l 


LEVEL UPPER 
320.2040 320.2040 
LEVEL UPPER 
~34.5557 +iNF 
0.0005 2.5000 
1.5248 2.0000 
0 NONOPT 


0 INFEASIBLE 

O UNBOUNDED 

0 ERRORS 
1.525 
5.000000E~4 


-34.556 


MARGINAL 
1.0000 
MARGINAL 


~53.5028 
EPS 


The optimal solution has the highest value for the maximum 


likelihood function. 


for the maximum likelihood function is 68.9552. 


parameters are: 
1. alpha (a) 
2. beta (b) 


3. theta (t) 


In this example, the highest value 


1.4218 


0.1623 


= ~0.1535 
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The optimal 





APPENDIX F 


Shown below is the CDF of the logistic function. 


if t#0 


if t=0 


For a given percentile (call it P), one sets G(x,;) = P , obtaining: 


a, 1f t#0 


. (x, 5a) 
1 +e 


After appropriate mathematical manipulation, one can take the 


natural logarithm and obtain the following equation: 
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Solving for x, one obtains the following expression: 


erle 


(2 - b (In (5-1) if t>0 


1 = 
= i c 1 
(a b (in ( i. 1))} if t<0 


Similar calculations were done to obtain :; 


(a - b (In (4-1))) 
xX, =e - if t=0 


2 


By substituting in values for a, b, t, and P the predicted 


commuting distance can be calculated. 
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APPENDIX G 
TABLE 14 


WESTERN UNITED STATES: A COMPARISON OF ACTUAL VERSUS LOGISTIC 
TRANSFORMATION COMMUTING DISTANCES 


eT iene: 
Percentile | Percentile | Percentile 
| act. | est. | act. | est. | act. | est. 
frucson, az | 582 9.7| 3.4| 27.2| 19.3] 62.8| 42.4 
‘concord, cA __—| 140 | 19.2| 18.5| 30.2| 29.6| 42.5. 

/23.7| 22.7| 31.7| 34.2| 60.4| 50.9 
El Monte, CA 21a 
Long Beach, CA__| 184 
Redding, CA 

3 











Dublin, CA 





46.3 


U 


San Pablo, CA 
Upland, CA 
Vallejo, CA 


0.4 | 49.6. 
.3| 19.4] 32.2] 30.3] 49.4. 
-13.e| 15.3| 22.0] 23.1] 28.8 | 34.8 
39.0 
13.2 
4 ; 
1.0 
<7 


7) 


oO | > Mire |S fe 
me fo am oe 
Ww {0 W To [hr {hd 


38.8 


Van Nuys, CA 


Aurora, CO 


44.3 
/14.4| 15.3| 45.8| 29.8| 63.9 | 60.1 
19.6| 20.4] 55.4 | 39.4 
30.2 


Aurora, CO 


Denver, CO 





38 
Hayden Lake, ID 13 
Great Falls, MT 
Helena, MT 


9 : 
7| 25.6| 21.6| 34.2| 30.1| 40.6 | 47.3 
| 6.0] 8.2] 19.2| 13.5] 57.6 

ja1.7| 14.6] 54.9] 2¢.6| 72.3 | 79.1 
9.4] 19.2] 33.9] 22.9] 46.0] 48-0 


~ 
Nh 
> 


Las Cruces, NM 188 
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TABLE 14 (CONTINUED) 


75th 90th 
Percentile | Percentile 


5.2 miles 8.6 miles 
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TABLE 15 


SOUTHERN UNITED STATES: A COMPARISON OF ACTUAL VERSUS LOGISTIC 
TRANSFORMATION COMMUTING DISTANCES 


50th 75th 90th 
Percentile | Percentile | Percentile 


est. 
/26.4| 29.1| 47.0| 42.1] 48.7 
-13.7| 13.6| 32.2] 26.5| 45.6| 48.7 
64.5 
1 BEE Wt) MEL oe 6744 
403 | 9.2| 9.4| 22.4| a7.3| 43.2] 35-5 
120 24.7 
365 
121 lo4.e] 23.1. : 
| 


Center Location 


Decatur, AL 
Dothan, AL 
East Camden, AR 





El Dorado, AR 
Little Rock, AR 
West Memphis, AR 
Seaford, DE 


VYallahassee, FL 


Bardstown, KY 
Slidell, LA 
Hattiesburg, MS 





56 
17 
87 
83 
5 32.2 4964 
63 38.6 51.6 
10 
0 
02 
55 
7 
52 


st. 
. 40.5. 
110 | 20.4 81.9 
Jackson, MS | 6.0] 7.4] 22.3] 23.6] 52.8 | 33-0 
et 


2 36.8 65.7 
236 18.3 34.2 


59.2 
/15.6| 29.7] 26.8| 44.6| 47.7 


: Greensboro, NC 


Wilmington, NC 
Dewey, OK 


3 /as.s| 18.4] 25.2| 28.2. 47.4 | 41.5 
152_| 6.4] 6.2] 9.2] 20.2] 37.2| 20.1 
1080 19.1 31.6 


1 
1 
9 
Greenwood, SC 
8 


Carrollton, TX 
Ft. Bliss, TX 
Houston, TX 
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TABLE 15 (CONTINUED) 


n 50th 75th 
percentile | percentile 

|__| act. | est. | act. | est. | act 

Pacis, TX 42.5 73.4 

[victoria, mx | 115 | 24.2 | 40.0| 58.8 

[Ft. story, va__| 681 | 15.5| 15.6| 21.2] 22.2] 31.8 | 32.9 

Richmond, VA 

Salem, VA 10.7 | 12.4 | 20.0 

Charleston, WV | S29 | 35.2] 13.7] 28.2| 24.0] 37.6] 41.1. 


Mean Absolute 1.8 miles 4.5 miles 7.5 miles 
Deviation = 









90th 
percentile 
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TABLE 16 


MIDWESTERN UNITED STATES: A COMPARISON OF ACTUAL VERSUS LOGISTIC 
TRANSFORMATION COMMUTING DISTANCES 
90th 


Center Location 50th 75th 
Percentile Percentile | Percentile 


pf acts | est. | act. | est. | act. 
james, TA 258 | 23.0] 23.7] 35.4] 29.6| 72.1] 65.5. 
‘council plutts, 1a |166| 10.0| 9.e| 20.4] 16.8| 38.4| 34.9. 
becorah, 18 | ag6 | 25.3| 27.4] 52.2] 45.3| 70.2] 70.9. 
Dubuque, TA 244] 8.2 10.0] 29.2| 20.2] 69.2] 54.3 
545 
/20.7| 22.2{ 24.8] 26.0] 28.6] 24.9 

3.5 

02 










| 24.0| 23.5| 42.9| 40.5] 67.4| 66.0 
/a3.2| 14.2] 34.4] 25.3 | 40.5] 40.2 
/ 6.2 7.0| 15.2] 12.4] 22.5] 23.3 
307 | 6.0] 7.3] 2.8] 22.0] 28.0] 27.5 


oie te 
N }Oo Ir 
mip te 


23.6| 20.5 
19.4 
|a9.9| 17a] 97.6] 34.0] 57.4] 64.4 
10.4| 9.4] 16.3] 18.0] 46.7] 36.2 
Fergus Falls, MN 122 


Bay City, MI 


Grand Rapids, MI 
Duluth, MN 






MM TH |W [Rk 
m 10 1O [WwW 
J jo Tn To 


columbia, Mo 256 | 5.5 


Washington, wo | aaa | 21.2] 21.0] 46.6] 35.7] 58.3] 60.1. 
Columbus, OH 713 8.5} 10.2} 35.8; 22.1 58.5] 53.6 





TABLE 16 (CONTINUED) 


od 
Percentile | Percentile {| Percentile 

__|act.| est. | act. | est. | act. | est. 
Toledo, OH 
Fond du tac, wr__| 160 | 4.0| 6-9] 20.6| 12.0| 28.1| 20.8 | 


Mean Absolute 1.3 miles 5.2 miles 


Deviation 
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NORTHEASTERN UNITED STATES: 


Center Location 








Huntingdon, PA 


Mean Absolute 
Deviation 





Percentile 


Lact. 
66.3 
17.9 
41.7 


act. 
7.5 
13.1 28.8 


st 
| 7s] 7.7] 16.5] 14.3 | 
5 


539 
135 
575 
883 
255 
112 
118 
474 
163 
152 15.5 
104 35.2 
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TABLE 17 


A COMPARISON OF ACTUAL VERSUS LOGISTIC 
TRANSFORMATION COMMUTING DISTANCES 


90th 
Percentile 


66.3, 


75th 
Percentile 


50th 


est. 


72.4 


/31.7| 36.2] 64.5 


30.7 
13.4 


42.7, 
28.8 


29.3 


st. 
| 4.o| 7. 
_28.0| 34.3| 36.8| 45.6| 46.6 
17.8| 25.1| 26.8] 38.2| 39.0 
| 20.9| 29.7] 32.2] 40.9] 47.2 
/15.5| 25.6| 24.1] 32.2] 35.9 
| 35.2| 36.3] 52.2 54.3 
| 20.4] 36.9 | 37.9 


41.0 
31.2 
21.4 
31.1 


1.1 miles 2.4 miles 3.6 miles 
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APPENDIX H 

A GAMS program was used to numerically solve the MLE's (n,s, 
and t) using the normal distribution for every reserve center. The 
input was formatted by using a fortran program. An example of the 
GAMS program is shown below. The reserve center depicted is 
located at Bakersfield, California. 

PROGRAM 

SOFFSYMLIST OFFSYMXREF 
OPTIONS SOLPRINT = On, LIMCOL = 0, LIMROW = 0; 


sets 


Sinclude file set; (this file assigns a number to each commuter; 
in this example, the numbers are from 1 to 149) 


scalar count 


Sinclude file count; (this file gives a scalar for the number of 
commuters to the reserve center; in this 
example, count equals 149) 


parameter X(T) 


Sinclude file out; (this file is a listing of the commute 
distances) 


variable 


value (this is the value of the maximum likelihood function) 
mu 

sigma 

th; (this is theta, the value of the power transformation) 


positive variable sigma; 


th.lo = 0.001; (these 5 lines establish bounds and set starting 
th.up ~ 1.000; values for theta, mu, and sigma; sigma can only 
mu.. 1 = 1.000; be a positve value; mu can take on any positve 
Sigma.l1 = 5.000; or negative value) 
th.l = 0.250; 
value.1l = 0.000; (initial value of the function is equal to zero) 
equation 

NORMAL; 
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NORMAL... value =e= 
((-0.5*count) *LOG(2*3.1459)) - (count) *LOG(sigma) 
- (0.5/(sigma**2) ) *(SUM(I, (ABS ((((X(I) **th) -1) /th) -mu) **2) ) ) 
+ (th-1)*SUM(I,LOG(X(I))); 


Model Stat /NORMAL/; 

Solve Stat using DNLP maximizing value; this directs the solver 
to find the maximum value for the function 
NORMAL using the Discontinuous Non-Linear 
Program algorithm) 


Display mu.1, Sigma.1, th.1, value.1; (this displays the values 
for mu, Sigma, theta, and the function NORMAL) 


th.lo = -0.950; (the value of theta is modified so that the 
theup = -0.001; function is evaluated from the left side 
mu. 1 = 1.000; of zero) 

Sigma.l1 = 5.000; 

th.l = ~0.200; 

value.l = 0.000; 


Solve Stat using DNLP maximizing value; 


Display mu.l, sigma.1, th.1, value.1; 


mu.l = 3.000; 
Sigma.l1 = 5.000; 
value.l = 0.000; 
equation 
NORMAL2 ; 


NORMAL2.. value =e= 
((-0.5*count) *LOG(2*3.1459)) - (count) *LOG(sigma) 
- ((0.5/(Ssigma**2) )*(SUM(I, (ABS (LOG(X(I) ) -mu) **2) ) ) ) 
- SUM(I,LOG(X(I))); 
Model Stat2 /NORMAL2/; 
Solve Stat2 using DNLP maximizing value; 


Display mu.1l, sigma.1, value.l1; (the value of theta equals zero) 
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LOWE 


EQU NORMAL 183.1801 


LOWER 
VAR Value ~INF 
VAR mu -INF 
VAR sigma : 
VAR th 0.0010 


REPORT SUMMARY: 


VARIABLE mu.l 
VARIABLE sigma.1] 
VARIABLE th.l 


VARIABLE value.l 


R 


Output 
LEVEL UPPER 


183.1801 183.1801 


LEVEL UPPER 
74.1694 +INF 
—-2.1459 +INF 

12519 +INF 

0.0010 1.0000 

0 NONOPT 


O INFEASIBLE 
0 UNBOUNDED 
0 ERRORS 
“2.146 

1.258 

0.001 


74.169 


MARGINAL 
1.0000 
MARGINAL 
EPS 
EPS 
~24.8407 





LOWER 


EQU NORMAL 183.1801 


LOWER 
VAR Value ~INF 
VAR mu ~INF 
VAR sigma : 
VAR th -0.9500 


REPORT SUMMARY: 


VARIABLE mu.l 
VARIABLE sigma.1 
VARIABLE th.1l 


VARIABLE value.l 


LEVEL UPPER 


183.1801 183.1801 


LEVEL UPPER 
75.9896 +INF 
-2.6778 +INF 

1.7013 +INF 
-0.1452 ~0.0010 

0 NONOPT 


QO INFEASIBLE 
O UNBOUNDED 
0 ERRORS 
-2.678 

1.701 
-0.145 
75.990 
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MARGINAL 
1.0000 
MARGINAL 
EPS 
EPS 
EPS 


LOWER LEVEL UPPER MARGINAL 


EQU NORMAL2 183.1801 183.1801 183.1801 1.0000 

LOWER LEVEL UPPER MARGINAL 
VAR value -INF 74.1942 +INF : 
VAR mu ~INF -2.1490 +INF EPS 
VAR sigma ; 1.2604 +INF EPS 
REPORT SUMMARY: 0 NONOPT 


O INFEASIBLE 
O UNBOUNDED 


0 ERRORS 
VARIABLE mu.1l = 2.149 
VARIABLE sigma.1 = 1.260 
VARIABLE value.l = 74.194 


The optimal solution has the highest value for the maximum 
likelihood function. In this example, the highest value 
for the maximum likelihood function is 75.9896. The optimal 


parameters are: 


1. mu (m) = =-2.6778 
2. Sigma (Ss) = 1.7013 
3. theta (t) = -0.1452 
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APPENDIX I 
TABLE 18 


SOUTHERN UNITED STATES: NORMAL DISTRIBUTION 


Center Location Le 
Decatur, AL -0.859 0 
Dothan, AL ~1.746 0.1455 
East Camden, AR -1.488 0.0746 
El Dorado, AR -1,556 0.0592 
Little Rock, AR -3.101 -0.2320 


t 
0.6430 


0.227 


0.790 


West Memphis, AR -4.457 1.791 

Seaford, DE “1.273 | 0.442 
Tallahassee, FL -2.378 -0.1000 
Bardstown, KY “1.226 0.2745 
Slidell, LA -0.814 0.7178 
Hattiesburg, MS -1.600 0.0353 
Jackson, MS -0.5512 
Greenskoro, NC -1.526 0.0461 
Wilmington, NC ~3.941 -0.4976 


Dewey, OK ~0.839 0.6968 
Greenwood, SC “2.210 0.1742 
Carrollton, TX -1.440 0.2036 
Ft. Bliss, TX 79906 ~6.036 -0.4968 
Houston, TX -2.397 | o.950 0.0984 
~1.264 0.2127 










Paris, TX 
Victoria, TX 77901 -1.397 0.774 0.1760 
Ft. Story, VA 23459 ~2.187 0.776 -0.1648 


~3.568 2.190 -0.3207 
-2.287 0.998 -0.1002 
~1.821 0.730 0.0849 


Richmond, VA 23220 
Salem, VA 24153 
Charleston, WV 25313 


© 
=) 


TABLE 19 


WESTERN UNITED STATES: NORMAL DISTRIBUTION 


ee ae: 
-3.306 1.804 -0.3259 
-2.678 1.701 -0.1452 


-2.003 1.144 ~0.1639 


“1.312 0.463 0.3310 
-1.404 | 0.590 | 0.0712 


“1.692 


“1.671 0.662 0.0879 


Center Location Zipcode 
85713 
93301 
95926 
94519 
94566 
91733 
90822 
96003 
94806 
91786 
94589 
91403 
80011 
80045 
80225 
83835 


Tucson, AZ 
Bakersfield, CA 
Chico, CA 
Concord, CA 
Dublin, CA 

El Monte, CA 
Long Beach, CA 
Redding, CA 
San Pablo, CA 
Upland, CA 
Vallejo, CA 


-0.0853 
-1.275 0.3238 
~1.721 0.1125 
-2.165 -0.0896 
-1.628 0.0980 
-2.682 -0.2316 
-1.954 -0.0560 
-2.663 -0.1760 


~2.250 -0.4708 


Van Nuys, CA 


Aurora, CO 
Aurora, CO 
Denver, CO 


Hayden Lake, ID 





Great Falls, MT 59403 ~3.290 ~0.2581 
Helena, MT 59601 -2.965 1.976 -0.4470 
Las Cruces, NM 88001 -3.328 1.833 -0.4724 


Reno, NV 


89502 
97402 
84321 
84407 


~6.144 3.278 ~0.5450 
Eugene, OR =3:i918 1.879 ~0.4468 


-5.838 3.198 -0.6331 





Logan, UT 


-1.980 0.923 ~0.0028 
-2.833 1.175 ~0.1651 
=Jei1/ 1.784 -0.3118 


Ogden, UT 


98404 
Spekane, WA 99216 
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TABLE 20 


MIDWESTERN UNITED STATES: NORMAL DISTRIBUTION 


Suen igasice _ugseie [as 
1.164 
[council Bluffs, ta_| 51501_| -6.e19 | 3.764 | -0.8476_ 
Granite city, 1 | 62080 | -4.016 | 1.662 _| -0.5120_ 
fquincy, 1) |__62902_| ~2.209 | 0.662 | 0.1910_ 
‘springfield, Ib . ) 




















| acsoo | -2.248 | 2.122 | -0.1505_ 


62703 “2.403 1.325 -0.1145 


-5.086 | 3.094 
{0.780 | 0.0290_ 
| -a.c6r_| 2.479 | -0.6222_ 

0.915 









Toledo, OH 
Beloit, WI 
Fond du Lac, WI 54935 
Wausau, WI 54401 
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TABLE 21 


NORTHEASTERN UNITED STATES: NORMAL DISTRIBUTION 


zipcode | om | is | 

0608s 

06457 

06120 
. 








02210 1.270 
01843 
02134 0.5463 
04404 
Dexter, ME 04930 
Pedrickton, NJ 08067 





14424 
13617 
12932 
12801 
15601 
Erie, PA 16504 
0.1199 
0.1572 
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APPENDIX J 


Shown below is the normal CDF. 


V4 -(vy-m) \2 
P= | : ‘ 7 ) av, 
“.V¥on S 





where P is the percentile commute distance. 
A statistical package called GRAFSTAT has a function labelled 
NORICDF. By entering the values for m,s, and the desired 
percentile (0.5 for the 50th percentile commute distance, 0.75 for 


the 75th percentile, etc.) into the function NORICDF, ore can 


obtain a value for v:. 


As discussed in Chapter III, the transformation selected is: 


x; co4 
ee 


In order to obtain the estimated commute distance, the 


transformation must be solved for x, . Therefore, the equation 


can be simplied to: 


i 
xX, = (tv;+1) 5 


The values of v, and t are substituted into the above equation to 


obtain the estimated commute distance. 
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APPENDIX K 
TABLE 22 


WESTERN UNITED STATES: A COMPARISON OF ACTUAL VERSUS NORMAL 
: TRANSFORMATION COMMUTING DISTANCES 


Center Location 50th 
Percentile 


est. 
| 9.7] 20.6] 27.2] 20.3, 
gu | 7.2 


n 










75th 
Percentile 


90th 
Percentile 





est. 


Tucson, AZ 


46.1 


582. 9 | 62.8 42.3 
Bakersfield, CA 7 61.8 
Chico, CA 
Concord, CA 140 | 19.2] a7.9| 30.2] 29.6] 42.5 
Dublin, CA 590 | 23 /22.8| 31.7| 35.2. 51.5 
El Monte, CA 639 27.9 
Long Beach, CA 184 40 

0 ; 


, Redding, CA 
San Pablo, CA 
Upland, CA 
Vallejo, CA 


| 25.2| 27.4 
| 10.4] 14.9] 34.2] 26.6, 


i7 
ae ; 
8 : 
2 
. 
5 
8 
.4 
8 ; 
37.4 | 39.2 


‘ek 


DN fk tre | PR Te 
Dm TR 1H fW IP 
O11 Fo. Te Fh 


Van Nuys, CA 


13.8 | 13.8 30.6 
38.7 
Aurora, CO : 44.2 
Aurora, CO 823 30.6 57.1 

7 


Ol 
\o 


Denver, CO 389 21.1 39.8 
40.6 

| 6.0] 9.2| 19.2] 17.5| 57.6 | 34.6 
67.4 
188 teu 13.5 33.9] 23.3] 46.0] 44.7 





Hayden Lake, ID 
Great Falls, MT 
Helena, MT 


bm fR 
N 7 UO 
> |] 


Las Cruces, NM 








TABLE 22 (CONTINUED) 
75th 


Center Location 50th 90th 
Percentile Percentile Percentile 


fact. | est. | act. | est. | act. | est. 
Reno, WV | 24a] 5.2| 6.7] 12.0] 12.2] 57.9] 26.6 
feugene, oR | 186 | 12.0] 12.1| 22.6] 22.1] 52.7| 40.6 
togan, ur 107 |_8.9| 8.7] 20.9] a5.0| 37.5] 32.0 
fogaen, vr | 403 | 20.0| 13.9] 35.5] 25.8] 42.4 45.2 
‘Tacoma, WA | 326| 8.2] 9.8| 20.3| 27.2] 32.7] 30.2 
[spokane, Wi. ___| 294 | 20.6] 12.3] 27.4] 22.3| 66.2] 47.8 


Mean Absolute 1.7 miles 4.7 miles 8.5 miles 
Deviation 
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TABLE 23 


SOUTHERN UNITED STATES: A COMPARISON OF ACTUAL VERSUS NORMAL 
TRANSFORMATION COMMUTING DISTANCES 
75th 


Center Location 50th 90th 
Percentile | Percentile | Percentile 


oo fact. | est. | act. 
fae ar ase fa a] 79.2 
117 
259 











-13.7[ 13.9] 32.2] 26.3| 45.6| 46.0 

22.5] 20.7| 38.3] 37.0] 62.7] 62.3 

‘El Dorado, aR | 187 | 26.6] 19.6| 42.7| 36.6 63.3 

Little Rock, aR__| 403 | 9.2] 9.7 

West Memphis, aR_| 120 | 12.1] 11.4{ 16.5| 16.4] 28.2 | 24.9 
7 





Bardstown, xv | 121 | 24.8| 22.4| 32.2| 34.7| 52.2 
31.5| 29.4| 38.6] 42.0] 49.2 | 52.3 
20.4 
[Jackson Ms___—|_509 | 6.0| 7.7| 22.3| 14.3] 52.8| 33.4 
‘Greensboro, nc___| 290 | 18.8| 20.6| 48.4| 38.2 
18.6| 44.4 | 33.7, 
Dewey, OK 102 

155 | 16.6| 15.4| 29.7| 26.5] 44.6 
: Houston, TX 12.4 3043) | 3255 





27 
ol 


Center Location 


Paris, TX 


Victoria, TX 
Ft. Story, VA 
Richmond, VA 


Salem, VA 
Charleston, WV 


Mean Absolute 
Deviation 





TABLE 23 (CONTINUED) 


50th 75th 
percentile | percentile 


_act.| est. | act. | est. 
6 


- | est. | act. | est. 
29.7| 22.9| 42.5| 42.3. 
EE Oe 
15.4 
9.3| 25.3| 20.2 
12.8 22.4 
| 28.2 | 


2.1 miles 


147 
115 
681 
573 
320 

29 


5 





percentile 
77.2 | 69.0. 
65.7 
33.7 
71.7] 50.6 


7.1 miles 


: 





, oe 


: Lake Station, IN 
i 


TABLE 24 


MIDWESTERN UNITED STATES: A COMPARISON OF ACTUAL VERSUS NORMAL 


TRANSFORMATION COMMUTING DISTANCES 
Center Location 







50th 75th 90th 
Percentile Percentile Percentile 


Lact. | est. | act.| est. | act. 
. 31.1 
[20.6 


ol 


m |e a 
m [0 
m [oa [Hn | o 


61.1 


act. 
| 
16 20.4 | 
_25.3| 32.0| 52.2 74.6 
| s.2| 10.3] 29.2 


69.1] 51.0 





Chicago, IL 


545 7.2 11.4 17.9 
Granite City, IL 11.3 16.9 | 28.6| 26.7 


63.1 


0 
5.3 50.4 
39.9, 


fuiney, 1b | 322 | 24.0] 22.8] 42.9] 39.9] 67.4 
58.7 


48.5 
21.5] 22.0 


11.8 18.2 


nN 
co 
~J 


32.7 | 24.8 
‘Fort Wayne, IN| 522 | 13.2| 24.4{ 34.4| 26.0 
[cary, IN| 208 


N 6 
Scottsburg, IN 23 40.8 
MI 9 
f 


2 


est 
16. 
6 td] 
ay City, 19.4 52.2 
Grand Rapids, ur__| 508 | 19.9 58.8 
Duluth, MN 2 9.3 . 
Fergus Falls, MN 


Columbia, MO 


10.4| 9.3] 16.3| 18.5] 46.7] 36.6 

33.9 | 51.0 62.9 
s.5| s.9| 2e.4| 17.1] 47.4| 42.8 
21.2] 21.2] 46.6| 36.5| 58.3| 59.3 
58.5 | 55.1 


306 
08 
47 

122 | 34.4 


Washington, MO 
Columbus, OH 
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TABLE 24 (CONTINUED) 


Eloi alin 
percentile percentile | percentile 
| act. | e act. | est. 
ffoledo, on | 495 | 6.5. 45.7 

0-4 


S 
7 
Fond du Lac, WI | 160 | aco 
225 | 18.3 


Mean Absolute 1.5 miles 
Deviation | 


ee eo ree a ae Te Re 


Center Location 


Lest. Lact. 
7.7 0] 45.7 | 39.4 
22.1 
| 6.0] 20.6] 9.9] 28.2] 21.2, 
56.6] 56.4 





TABLE 25 


NORTHEASTERN UNITED STATES: A COMPARISON OF ACTUAL VERSUS NORMAL 


TRANSFORMATION COMMUTING DISTANCES 
Center Location Pe 
fact | est. | act. | est. | act. | est. 
East Windsor, CT 53.3 
Middletown, cr__| 135 | 18.4| 18.1| 312.7] 36.4| 94.5 | 67.7 
West Hartford, CT 41.9 
9| 4. 6. 
7 


50th 75th 90th 
Percentile Percentile Percentile 


act 
7 


9 
Bangor, mE | 2 | 15.0] 15.4 
Dexter, ME ais | 

474 
4 


Pedrickton, NJ 


15.0| 15.4] 34.6| 28.4] 61.0| 55.7 


40.1 


Roslindale, MA . 9 13.6 


112 
118 : : 
474 
‘canadaigua, ny | 163 | 25.2] 20.1| 29.7] 32.1] 40.9 
16.9 S22) hb 35.7 
23.7 35.2 55.5 
29.0, 


47.1 


18.7| 20.3| 28.8 | 31.6 
7 : 
_29.3| 27.6] 34.3| 37.8] 45.6] 49.2 


36.5 


3 | 26.0| 20.5| 36.9 
30.4 
6.0 38.6 
32.9 
22.7 
funiontown, pa | 350 | 13.4] 22.3| 20.6] 19.3] 30.6] 30.3 


Mean Absolute 1.2 miles 2.6 miles 4.4 miles 
Deviation 







22 


Glens Falls, NY 
1 
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APPENDIX L 
Listed below are the nine census regions that are used in the 
multiple regression model. Also shown are the states that comprise 


each region. 


Region States 
1. Pacific Washington, Oregon, California 
2. Mountain Montana, Idaho, Nevada, Utah, 


Arizona, Colorado, Wyoming, 
New Mexico 


3. West North Central North Dakota, South Dakota, Iowa, 
Nebraska, Kansas, Minnesota, 
Missouri 


4. West South Central Texas, Oklahoma, Arkansas, 
Louisiana 


5. East North Central Wisconsin, Illinois, Michigan, 
Indiana, Ohio 


6. East South Centra) Kentucky, Tennessee, Mississippi, 
Alabama 


7. New England Maine, Vermont, Massachusetts, 
Connecticut, Rhode Island, 
New Hampshire 
8. Mid-Atlantic New York, Pennslyvania, New Jersey 
9. South Atlantic Maryland, West Virginia, Virginia, 


North Carolina, South Carolina, 
Georgia, Florida, Delaware 
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APPENDIX M 
Listed below are the types of units located at the reserve 
centers and the category to which they are assigned. 
Category Unit Type 


1. Combat Airborne, Armour, Field Artillery, 
Infantry, Special Forces 


2. Combat Support Headquarters, Air Defense, Aviation, 
Chemical, Engineers, Logistical 
Command, Military Intelligence, 
Military Police, Signal 


3. Combat Service Adjutant General, Civil Affairs, 
Support Finance, Judge Advocate General, 
Ordinance, Psychological Operations, 
Quartermaster, Transportation, 
Training, Schools 


4. Medical Medical 
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APPENDIX N 
TABLE 26 


SUMMARY STATISTICS FROM THE MULTIPLE 
REGRESSION ANALYSIS OF VARIANCE TABLES 


ee ee ee 
0.1614 
[Error Sum of Squares | 362.2523 | 90.6151 | 7.9743 
Total Sum of Squares 93.3035 
Model Mean Square 0.2569 

| 2.7400 | 


1.7400 | 2.2000 | 2.7000 
0.1223 | 0.0505 | 0.0189 





a — ) 








APPENDIX O 
TABLE 27 


PARAMETER ESTIMATES FROM THE 
MULTIPLE REGRESSION MODEL FOR T 


Variable Parameter Standard 
Estimate Error of 
Estimate 


Intercept Term -0.7518174 0.3471955 
Male III-A Fopulation 9.89541E-07 0.0000006 0.0916 


Percentage of Males at 0.6258805 0.2877690 
the Reserve Center 

Civilian 0.0047917 0.0165479 
Wages/Military Wages | 


Unemployment 0.0123455 0.0178909 


Percentage of Prior -0.3333634 0.3445064 
Service Personnel __ 


Propensity to Commute 0.0063553 0.0024008 0.0097 
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TABLE 28 


PARAMETER ESTIMATES FROM THE 
MULTIPLE REGRESSION MODEL FOR M 


Variable Parameter Standard 
Estimate Error of 
Estimate + 


Intercept Term -~2.7087299 2.3393631 
Male III-A Population 4.48160E-06 0.0000039 


Percentage of Males at 3.6032658 1.9389543 
the Reserve Center 

Civilian ~0.2350680 0.1114978 
Wages/Military Wages 


Unemployment 


Percentage of Prior 
Service Personnel 


Propensity to Commute 0.0291351 0.0161764 0.0753 
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TABLE 29 


PARAMETER ESTIMATES FROM THE 
MULTIPLE REGRESSION MODEL FOR S 


Variable Parameter Standard Prob. > 
Estimate Error of | t | 





Estimate 


Intercept Term 3.2768767 1.1035709 0.0039 
Male III-A Population -3.96858E-06 0.0000018 0.0342 


Percentage of Males at -1.5044185 0.9146821 0.1038 
the Reserve Center 

Civilian 0.0213407 0.0525979 0.6860 
Wages/Military Wages 


Unemployment 0.0644819 0.0568669 0.2601 


Percentage of Prior -0.5948843 1.0950236 0.5884 
Service Personnel 


Propensity to Commute -0.0184917 0.0076310 0.0175 


APPENDIX P 
TABLE 30 


PARAMETER SUM OF SQUARES FROM THE MULTIPLE REGRESSION 
MODEL FOR S INCLUDING THE NINE CENSUS REGIONS 









Prob. > 
[t| 


Standara 
Error of 
Estimate 


Intercept Term 3.6295838 0.0000021 0.0962 
Male III-A Population ~3.54941E-06 1.2976608 0.0065 


Percentage of Males at ~1.8115175 0.9680104 0.0651 
the Reserve Center 

Civilian -0.0104021 0.0583026 0.8589 
Wages/Military Wages 


Unemployment -0.0698944 0.0721239 0.3356 


Percentage of Prior -0.8899148 1.4312514 - 0.5360 
Service Personnel 


Propensity to Commute 0.0081646 0.0431 
Census Regions 
0.6406002 0.1404 
0.1964427 0.7169 
0.0473899 0.9191 
0.8647664 0.0532 


0.1205138 0.5495573 0.8270 
0.2712929 0.5105324 0.5967 


South Atlantic -0.0913791 0.5097881 0.8582 
West North Central 0.5724342 0.5009444 0.2567 
West South Central 0.0000000 Base Case 


Overall Analysis of Variance (ANOVA) statistics are shown below. 


Parameter 
Estimate 


Variable 








East North Central 
East South Central 
Mid-Atlantic 


Mountain 






New England 


Pacific 


R-Square = 0.2285 
Model S.S.. = 21.3213 
Error S.S. = 71.9823 
F Value = 1.61 
Prob. > F = 0.0962 


TABLE 31 


PARAMETER SUM OF SQUARES FROM THE MULITPLE REGRESSION 
MODEL FOR M INCLUDING THE NINE CENSUS REGIONS 


v Variable Parameter Standard Prob. > 
Estimate Error of lt | 
Estimate 
: Intercept Term 0.9144551 2.6660525 0.7325 
Male III-A Population 9.85805E-07 0.0000043 0.8205 


3.9223890 1.9887836 0.0522 
Civilian -0.3825306 0.1197831 0.0020 
Wages/Military Wages 


Unemployment -0.1552018 0.1481790 0.2982 


Service Personnel 

Propensity to Commute 0.0425049 
Census Regions 
Bast North Central 













Percentage of Males at 
the Reserve Center 



























0.3884 
0.1774 


Overall Analysis of Variance (ANOVA) statistics are shown below. 


R-Square = 0.2538 
Model S.S. = 103.3648 
g Error &.S. = 303.8369 
F Value = 1.85 
Prob. > F = 0.0466 


TABLE 32 


PARAMETER SUM OF SQUARES FROM THE MULTIPLE REGRESSION 
MODEL FOR T INCLUDING THE NINE CENSUS REGIONS 


Variable Parameter Standard Prob. > 
Estimate Error of lt | 
Estimate 


Intercept Term -0.5454899 0.3858434 
Male III-A Population 6.07259E~06 0.0000006 


Percentage of Males at 0.7014002 0.2878259 
the Reserve Center 


Civilian 0.0054210 0.0173356 
Wages/Military Wages 


Unemployment 0.0003004 0.0214452 


Percentage of Prior -0.4388938 0.4255649 
Service Personnel 


Propensity to Commute 0.0072217 0.0024277 





0.1615 
0.3365 
0.0172 


0.7554 


0.9889 
0.3057 


0.0039 








East North Central 


East South Central 


Census Regions 


-0.2368959 0.1278481 0.0678 
-0.0296627 0.1605147 0.8539 






Mid-Atlantic -0.1412891 0.1382125 0.3099 
-0.4817537_| 0.1309262 | 0.0004 


New England 





Pacific 


South Atlantic ~0.1563604 0.1515792 0.3056 


0.3179 
0.6416 


~0.1642686 0.1634041 
~0.0709449 0.1518005 





West North Central ~0.2495216 0.1489496 0.0980 
West South Central 0.0000000 Base Case Base Case 


Overall Analysis 
R-Square 
Model S.S. 
Error S.S. 
F Value 


Prob. > F 


of Variance (ANOVA) statistics are shown below. 


0.3316 


3.1570 


il 


6.3639 


= 2.69 


0.0030 
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TABLE 33 


PARAMETER SUM OF SQUARES FROM THE MULTIPLE REGRESSION 


Variable 


Intercept Term 


Civilian 


MODEL FOR T INCLUDING UNIT TYPE 


Parameter Standard Prob. > 
Estimate Error of 
Estimate 


-~0.7251097 0.3820098 0.0612 





6) 


Percentage of Males at 0.4748144 0.4439196 0.2880 
the Reserve Center 


0.0065104 0.0173431 0.7084 
Wages/Military Wages 


Unemployment 


Service Personnel 


Percentage of Prior ~0.2009020 0.4044274 


Male III-A Population 9.09502E-07 0.0000006 0.1321 


0.0093589 0.0186847 


0.6178 
0.6207 





Propensity to Commute 0.0061539 0.0024986 0.0159 


Combat 
Combat Support 


Combat Service Support 0.0387408 0.1327390 0.7711 


Medical 


Overall Analysis 
R-Square 
Model S.S. 
Error S.S. 
F Value 


Prob. > F 


Unit Type 


0.0906429 0.1632313 


0.5802 


0.1067692 0.1314565 


0.4191 


0.0000000 Base Case 


of Variance (ANOVA) statistics are shown below. 


0.1704 


1.6224 


= 7.8986 


1.85 


0.0719 


101 





TABLE 34 


PARAMETER SUM OF SQUARES FROM THE MULTIPLE REGRESSION 
MODEL FOR M INCLUDING UNIT TYPE 




















Standard 
Error of 
Estimate 


Intercept Term -2.9060647 2.5470613 0.2573 
Male III-A Population 4.93705E-06 0.0000040 0.2191 


Percentage of Males at 4.4396847 2.9598470 0.1375 ~ 
the Reserve Center . 

Civilian -0.23008624 0.1156354 0.0500 
Wages/Military Wages 


Unemployment 0.0108268 0.1245809 0.9310 


Service Personnel 
Unit Type 


Overall Analysis of Variance (ANOVA) statistics are shown below. 


Parameter 
Estimate 


Variable 








R-Square = 0.13769 
Model S.S. = 56.0662 
Error S.S. = 351.1355 
F Value = 1.44 
Prob. > F = 0.1863 
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TABLE 35 


PARAMETER SUM OF SQUARES FROM THE MULTIPLE REGRESSION 


MODEL FOR S INCLUDING UNIT TYPE 
Prob. > 
|t| 















Parameter 
Estimate 





Variable Standard 
Estimate 


Error of 
Intercept Term 2.9771207 1.2133411 0.0163 
Male III-A Population -3.75635E-06 0.0000019 0.0513 


Percentage of Males at -0.8058185 1.4099795 0.5692 
the Reserve Center 

Civilian 0.0175577 0.0550851 0.7507 
Wages/Military Wages 

Unemployment ~0.0508008 0.0593465 0.3945 


Percentage of Prior -1.1647825 1.2845440 0.3672 


Service Personnel 
00079362 
Unit Type 
0.6522 


Overall Analysis of Variance (ANOVA) statistics are shown below. 



















Propensity to Conmute 





R-Square = 0.1459 
Model S.S. = 13.6213 
Error S.S. = 79.6823 
F Value = 1.54 

Prob. > F = 0.1486 
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